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ABSTRACT 

We propose to identify pulsar-wind bubbles (PWBs) as the environment in 
which the afterglow emission in at least some gamma-ray burst (GRB) sources 
originates. Such bubbles could naturally account for both the high fraction of 
the internal energy residing in relativistic electrons and positrons (e e ) and the 
high magnetic-to-internal energy ratio (e#) that have been inferred in a number 
of sources from an interpretation of the afterglow emission as synchrotron 
radiation. GRBs might occur within PWBs under a number of scenarios: in 
particular, in the supranova model of GRB formation a prolonged (months to 
years) period of intense pulsar-type wind from the GRB progenitor precedes the 
burst. Focusing on this scenario, we construct a simple model of the early-time 
structure of a plerionic supernova remnant (SNR), guided by recent results 
on the Crab and Vela SNRs. The model is based on the assumption of an 
“equipartition” upper bound on the electromagnetic-to-thermal pressure ratio 
in the bubble and takes into account synchrotron-radiation cooling. We argue 
that the effective upstream hydrogen number density for a relativistic shock 
propagating into the bubble is given by nH,equiv = [4p + ( B' + £') 2 / A-k\ / m p c 2 , 
where B' and S' are, respectively, the comoving magnetic and electric fields 
and p is the particle pressure. We show that, for plausible parameter values, 
^H,equiv spans the range inferred from spectral fits to GRB afterglows and 
that its radial profile varies within the bubble and may resemble a uniform 
interstellar medium or a stellar wind. We consider how the standard expressions 
for the characteristic synchrotron spectral quantities are modified when the 
afterglow-emitting shock propagates inside a PWB instead of in a uniform 
interstellar medium and demonstrate that the predictions for the empirically 
inferred values of e e and cb are compatible with the observations. Finally, we 
outline a self-consistent interpretation of the X-ray emission features detected in 
sources like GRB 991216 in the context of the supranova/PWB picture. 

Subject headings: gamma rays: bursts — MHD — pulsars: general — pulsars: 
individual (Crab Nebula, Vela Pulsar) — shock waves — supernova remnants 
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1. Introduction 


Gamma-ray burst (GRB) sources are commonly interpreted in terms of nonthermally 
emitting shocks associated with relativistic (and possibly highly collimated) outflows from 
stellar-mass black holes or strongly magnetized and rapidly rotating neutron stars (see, e.g., 
Piran 1999 and Meszaros 2001 for reviews). The prompt high-energy emission is thought 
to originate in the outflow itself, with the y-rays attributed to internal shocks within the 
flow and with the associated optical “flash” and radio “flare” emission ascribed to the 
reverse shock that is driven into the outflowing material as it starts to be decelerated by the 
inertia of the swept-up ambient gas. By contrast, the longer-term, lower-energy afterglow 
emission (see, e.g., |van Paradijs, Kouvcliotou, & Wijcrs 2000| for a review) is attributed to 
the forward shock that propagates into the ambient medium. The ambient gas is usually 
taken to be either the interstellar medium (ISM) of the host galaxy or a stellar wind from 
the GRB progenitor star. 

It appears that most of the observed emission from GRBs and their afterglows 
represents synchrotron radiation (e.g., [Sari, Piran, fe Narayan 1998| ; Panaitescu fe Meszaros 
19981 ; ISari fc Piran 19991 ; IGhcvaficr fc hi 200q ; |Lloyd fc Petrosian 20001) . In view of 


source-energetics considerations, the emission efficiency must be high. This implies that the 
ratio e e of the internal energy in relativistic electrons and positrons to the total internal 
energy density in the emission region is not much smaller than 1, and that the ratio e# 
of the magnetic-to-internal energy densities is not much smaller than e e . If the shocked 
gas consists of protons and electrons, then only moderately high (<; 0.1) values of e e may 
be expected even under optimal circumstances (e.g., |Bykov fc Meszaros 1996 ; [Potil fc 
Schlickciscr 2000|) . For e e to approach 1, it is probably necessary for the preshock gas 
to be composed primarily of e* pairs. A pair-dominated outflow is, in fact, a feature of 
certain GRB models (e.g., |Usov 1994; [Meszaros fe Rees 1997; |Grimsrud fe Wasserman 


1998|). Furthermore, the radiative efficiency of the reverse shock (and possibly also of the 


forward shock during the early afterglow phase) could be enhanced through pair creation 
by the high-energy photons comprising the gamma-ray pulse (e.g., IPhompson fe IVladau 
|2000| ; |Meszaros, Ramirez-Ruiz, fe Rees 200ip . There is, however, no natural way to account 
for large values of e e during the later phases of afterglows in a typical ISM or stellar-wind 
environment. 

It is in principle also possible to account for comparatively large values of in internal 
and reverse shocks by appealing to shock compression of magnetized outflows (e.g., |Spruit 


Paignc, fc Brcnkhahn 20UT| ; |Granot. fc Konigl 200 1| , hereafter GK). However, in the case of 
afterglows in the standard scenario, the highest values of €b that might be attained in this 
fashion (e.g., in a shock propagating into a magnetized wind from a progenitor star; see 
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Biernmnn fe Gassinelli 1993[) could at best account only for the low end of the actual range 
inferred in GRB afterglows (e# ,> 10 -5 ; e.g., |Panaitescu fc Kumar 200^ ). For example, 


one could not explain in this way the estimate e# ~ 0.01 — 0.1 (derived by model fitting 
of one of the most comprehensive spectral data sets obtained to date) in the GRB 970508 
afterglow (e.g., IWijcrs fc Galama 1999| ; IGranot, Piran, fc Sari 19991 ; jGhcvalicr fc Li 20001) fl 


As an alternative to compressional amplification of a preshock field, various proposals 
have been advanced for generating strong magnetic fields in the shocks themselves, but it 
is still unclear whether any of them can naturally account for a source like GRB 970508. 
For example, Medvedev & Loeb (1999) suggested that a two-stream instability behind the 
shock can generate fields that fluctuate on the very short scale of the plasma skin depth. 
However, the most likely value of e# predicted by this scheme is rather low (<^ 0.01), as 
is also the value of e e (e.g., [Pohl fe Schlickeiser 2000|); furthermore, questions have been 


raised about whether the fields will not, in fact, be damped on a similar microscopic scale 
(Gruzinov 1999). Thompson & Madau (2000) suggested that acceleration of the preshock 
gas by the prompt gamma-ray pulse photons would induce shearing motions that could 
significantly amplify the ambient magnetic field, ft is, however, unlikely that the preshock 
optical depth would be large enough for this effect to play a role for the comparatively 
large spatial scales (^ 10 1 2 ' cm) and low preshock densities (~ 0.03 — 3 cm -3 ) inferred 
for the GRB 970508 afterglow (see Meszaros et al. 2001). Blandford (2000), arguing by 
analogy with supernovae like Gas A, proposed that the afterglow emission in a source like 
GRB 970508 arises near the contact discontinuity that separates the swept-up ambient gas 
from the outflowing matter, where these two components can mix and interact. The large 
inferred magnetic field presumably originates in the central source and undergoes additional 
amplification in the turbulent interaction zone, but a quantitative model of this scenario 
has not yet been presented.^ 

In this paper we propose that the large values of eg and e e inferred in afterglows 
like GRB 970508 arise naturally if the outflow that gives rise to the gamma-ray pulse 


1 The inferred value of e e in this source is also fairly high (~ 0.1 — 0.6). 

2 A related idea was discussed by Smolsky & Usov (2000), who considered a magnetized, pulsar-type wind 
and suggested that it does not initially form a forward shock but rather that the oscillating currents in 
the wind front excite large-amplitude electromagnetic waves in the ambient medium and that high-energy 
electrons accelerated in the front radiate in the field of these waves. However, these authors still attributed 
afterglow emission detected more than a day after the ganrma-ray pulse (as was the case for GRB 970508) 
to a conventional forward shock that develops ahead of the wind front by that time. Additional ideas on 
how large-amplitude electromagnetic waves in Poynting flux-dominated outflows could lead to large values 
of £b in afterglows were outlined by Lyutikov & Blackman (2001). 
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expands into a pulsar-wind bubble (PWB). Such a bubble forms when the relativistic wind 
(consisting of relativistic particles and magnetic fields) that emanates from a pulsar shocks 
against the ambient gas and creates a “pulsar nebula,” whose structure is analogous to 
that of a stellar wind-inflated “interstellar bubble.” When a bubble of this type expands 
inside a supernova remnant (SNR), it gives rise to a “plerionic” SNR, of which the Crab 
and Vela remnants are prime examples (see, e.g., Chevalier 1998 for a review). GRBs can 
arise inside PWBs under a number of plausible scenarios, some of which have already been 
considered in the literature. For example, Gallant & Achterberg (1999) suggested that, if 
GRB outflows are formed in neutron-star binary mergers and expand into PWBs created 
by the progenitor pulsars, then acceleration of relativistic ions in the nebula by the forward 
shock could in principle account for the observed population of ultra-high-energy cosmic 
rays (UHECRs).]] Vietri & Stella (1998, 1999) presented a scenario for the origin of GRBs 
in which a rotationally supported “supramassive” neutron star (SMNS) forms either by a 
supernova explosion that is triggered by the collapse of a massive star or as a result of mass 
transfer in a low-mass X-ray binary (LMXB). In this picture, the neutron star loses angular 
momentum (and associated centrifugal support) through a pulsar-type wind until (on a 
time scale of several months to several years) it becomes unstable to collapse to a black hole 
(a process that, in turn, induces a GRB outflow). Vietri & Stella (1998, hereafter VS) noted 
the analogy between the proposed “supranova” remnants and plerionic SNRs, but they did 
not explicitly address the structure of SMNS wind nebulae and their implications to GRB 
afterglows.]] The afterglow sources observed to date are associated with “long” bursts (of 
duration ^ 2 s) and are often found within the projected image of the host galaxy. Such 
sources could plausibly arise in the collapse of (or the merger of a compact object with) a 
massive star (e.g., Woosley 2000), although an LMXB progenitor may also be consistent 
with the data ([Vietri fc Stella 19991) • In view of the evidence that at least some afterglow 
sources are located along the line of sight to a star-forming region (e.g., Meszaros 2001), 


3 The association of UHECRs with GRBs was originally proposed by Waxman (1995), Milgrom & Usov 
(1995), and Vietri (1995). Some difficulties with the simplest formulation of this idea were recently discussed 
by Stecker (2000) and Farrar & Piran (2000). UHECRs may, however, originate in winds from young, rapidly 
spinning and strongly magnetized neutron stars even if the latter are not linked to GRBs (see [Blasi, Epstein, 
fe Olinto 2000|). 


4 In a recent paper, Inouc, Guetta, fc Pacini 2001 discussed some observational consequences of the 
precursor plerion in the supranova scenario as well as the prospects for its direct detection. Their picture 
differs from ours in that they assume that the PWB disrupts and penetrates through the supernova ejecta 
shell, attaining a size that is about an order of magnitude larger than the SNR radius, and they associate the 
afterglow-emitting gas with entrained SNR fragments or the ambient ISM. By contrast, in our picture the 
PWB remains largely confined within the SNR shell, and the afterglow emission arises in the high-effective- 
density interior gas. 










5 


we adopt the supranova version of the SMNS scenario (VS) as the underlying framework 
for our discussion. 5 In §2 we estimate the physical parameters of SMNS winds and of 
supranova remnants in light of recent work on plerions, and we then model the structure 
of the resulting PWBs. In §3 we consider the expected properties of GRB afterglows that 
originate in such an environment. Our conclusions are summarized in §4. 


2. Pulsar-Wind Bubbles in Young Supernova Remnants 
2.1. The Supranova Scenario 


Supramassive neutron stars are general-relativistic equilibrium configurations of rapidly 
rotating neutron stars whose masses exceed the maximum mass of a nonrotating neutron 
star (e.g., |Cook, Shapiro, fc Tcukolsky 1994| ; |Salgado ct ah 1994p . A uniformly rotating 
SMNS that loses energy and angular momentum adiabatically while conserving its total 
baryon mass follows an evolutionary sequence that brings it to a point where it becomes 
unstable to axisymmetric perturbations, whereupon it undergoes a catastrophic collapse 
to a black hole. In their supranova model, VS postulated that the SMNS, which forms in 
the course of a supernova explosion of a massive star, is magnetized and loses energy and 
angular momentum through a pulsar-type wind|] The rate of energy loss can be estimated 
from the magnetic dipole-radiation formula 


Lir 




6 c 3 




ergs s 1 , 


( 1 ) 


where R* is the polar surface magnetic held, R * is the circumferential radius (neglecting 
the distinction between its equatorial and polar values in this approximation), and fl* is 
the (uniform) angular velocity (whose maximum value is ~ 2 x 10 4 s _1 ; e.g., |Haensel, 


'[ It has not yet been explicitly demonstrated that the supranova scenario can account for long bursts; in 
fact, it has even been suggested (Bottcher & Fryer 2001) that this model is most likely to produce short 
bursts. We note, however, that long bursts could in principle be generated in the course of the collapse 
of the SMNS to a black hole (see Kluzniak fe Ruderman 199§| ). Alternatively, if (as suggested by VS) the 
GRB outflow is produced after the collapse in a magnetized debris disk formed by the outer layers of the 


SMNS, then a long duration could be a consequence of a comparatively low disk viscosity (see, e.g., Popharn 


Woosley, k, Fryer 1999| and Ruffert fe Janka 1999| ) or of a magnetically mediated spin-up torque exerted by 
the black hole (van Putten & Ostriker 2001). 


6 As noted by VS, gravitational radiation, possibly associated with the excitation of r modes in the SMNS 
is an alternative loss mechanism. However, in view of the uncertainties involved in 


(e.g., |Andersson 1 

quantifying this process, we follow VS and neglect it in the ensuing discussion. 
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Lasota, & Zdunik 1999). [] The magnetic field amplitude in this estimate is normalized by 
the typical radio-pulsar value. This situation is to be distinguished from scenarios in which 
a magnetized rotator with a much stronger field (;> 10 15 G) is invoked to account for the 
GRB outflow itself (e.g., |Usov 19fTT| ; [Thompson 1994| ; |Blackman &; Yi 1998| ; |KluzniaOZ 


Ruderman 1998| ; ISpruit 1999| ; [Ruderman, Tao, fe Kluzniak 2QQ0Q . The initial neutron-star 
magnetic field might be amplified to such strengths by differential rotation (e.g., [Kluzniak 
|fc Ruderman 1998PI ) or through dynamo action (e.g., [Thompson fe Duncan 1993H ) . For the 
comparatively low field amplitudes adopted in the supranova scenario, the dynamical effect 
of the magnetic field on the SMNS structure should be negligible (see, e.g., [Bocquet et al 


19951) . 


The wind power L w consists of electromagnetic and particle contributions. The 
magnetic field is expected to have a dominant toroidal component, which scales with 
distance r from the center as 1/r. Under ideal-MHD conditions, the Poynting-to-particle 
energy flux ratio in the wind is given by 

^ > ( 2 ) 
47T p w c 2 

where B w is the field amplitude and p w is the rest-mass density (both measured in the fluid 
frame), and it remains constant along the flow after the terminal speed is reached. There 
has been a great deal of debate in the literature about the value of a w in relativistic pulsar 
outflows and about whether an ideal-MHD description is appropriate (see, e.g., Arons 1998 
and Chinch, Li, fc Bcgclman 1998| for discussions of this topic). For example, dynamical 
and radiative models of the Crab pulsar nebula have yielded preshock values of a w in the 
range ~ 2-5 x ICC 3 (e.g., [Arons 2002|) .| rfT | On the other hand, recent X-ray observations 


7 The magnetic dipole luminosity also scales as sin 2 (9qb, where 9qb is the angle between the rotation and 
dipole axes. However, as the spin-down torque of radio pulsars is evidently largely independent of the value 
of 9qb (e.g., Bhattacharya fe Srinivasan 1995| ), we did not include this factor in equation (|lj). 


®These authors also mention the possibility, suggested independently by VS in their outline of the 
supranova scenario, that differential rotation leading to a very strong field and possibly a GRB outflow 
could be induced in an SMNS when it starts to collapse after losing centrifugal support. 


9 These authors associate the dynamo action with neutrino flux-driven convection, which should occur if 
the stellar spin period 27t/H is shorter than the convective overturn time 10 3 F 39 1 ^ 3 s at the base of the 
convection zone, where F 39 is the neutrino heat flux in units of 10 39 ergs cm -2 s -1 , the expected value for a 
typical supernova explosion. Since the energy released in a supranova explosion may be significantly larger 
than in a typical supernova, and since the neutron-star mass is also higher in this case, the neutrino flux 
could be similarly larger and the convection condition might not be satisfied, which would be consistent with 
the assumption that the SMNS magnetic field does not exceed typical pulsar values. 


10 Begelman (1998), however, argued that a key underlying assumption of some of the dynamical estimates 
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of the Vela pulsar nebula have been interpreted as implying a w ~ 1 QHclfand, Gotthclf, fc 
Halpcrn 2001| ).p| In view of this uncertainty, and in order to examine the dependence of 
our model results on this parameter, we derive solutions for a w in the range 10 -3 — 1 . 

We also need to specify the Lorentz factor 7 ^ and the composition of the outflow. 
Spectral (e.g., [Kennel fe Coroniti 1984|) and optical brightness-distribution QGallant fc Arons 
|1994|) models of the Crab nebula have implied a current value of ~ 3 x 10 6 for 7 ^ upstream 
of the shock. However, although these fits account for the optical-through-gamma-ray 
observations, they do not explain the measured radio spectrum. I 11 a recent model, Atoyan 
(1999) interpreted the latter as being produced by a relic population of relativistic electrons 
that had been accelerated during the early years of the pulsar and that have subsequently 
lost most of their energy by radiation and adiabatic-expansion losses. Based on this 
interpretation, he argued that the Crab pulsar was born with a period of ~ 3 — 5 ms (as 
compared with previous estimates of ~ 19 ms)[] and initially had 7 ™ < 10 4 . In light of 
this result, we adopt 7 ™ = 10 4 as a fiducial value in our calculations: we assume that its 
magnitude is roughly the same in all objects and that it does not change significantly over 
the SMNS spin-down time. The pulsar outflow could consist of pairs as well as ions. In 
fact, by modeling the wind termination shock in the Crab nebula, Gallant & Arons (1994) 
inferred that the energy flux in ions is approximately twice that in pairs in that source, 
and we already mentioned (see §1) the suggestion by Gallant & Achterberg (1999) that 
UHECRs might be identified with heavy ions in GRB-associated PWBs. Nevertheless, for 
simplicity, we assume in our model that the SMNS wind is composed purely of e ± pairs.0 
In this case the wind power can be written as 


L w = 4vr(l + cr w )r 2 n w (r)'flP w 7n e c 3 , 


(3) 


where n w {r ) is the fluid-frame wind particle density at a radius r, m e is the electron mass, 


- that the magnetic field inside the shocked-wind bubble maintains a coherent, large-scale, toroidal structure 
— may not be valid, and he suggested that a w could be higher in this source. 

11 This interpretation was, however, questioned by Arons (2002), who suggested that a w in the Vela nebula 
is, instead, < 0.05; an even lower upper limit was given by Chevalier (2000), who proposed that a w < 10~ 4 
in this source. 


12 Independent arguments for why radio pulsars like the Crab and Vela were likely born with rotation 
periods <( 1 ms were recently given by Lai, Chernoff, & Cordes (2001) in the context of an interpretation of 
the apparent alignment of the spin axes, proper motion directions, and polarization vectors of the Crab and 
Vela pulsars. 


13 If the ion and pair energy fluxs are indeed comparable and these two components do not exchange energy 
efficiently, then the resulting bubbles would be approximated by the weakly cooling PWB solutions presented 
in §ED. 
















and (3 w c is the wind speed.Q 


The spin-down time of a rapidly rotating SMNS can be estimated as 

_A E mt fa \( M* \ 2 / R* \“ 6 / n* \~ 3 f B * \“ 2 
sd ~ 6 V0.5J ( 2.5 M 0 J V15knJ UO 4 s" 1 ] UO 12 G/ yr 


( 4 ) 


(see VS), where A E rot = aGMfQ*/2c is the portion of the rotational energy of an SMNS 
of mass M* and angular velocity hi* that needs to be lost before it becomes unstable to 
collapse.^] The basic time scale is determined by the underlying physical picture of a 
magnetized neutron star in which a significant fraction of the binding energy is invested 
in rotation (which is uniform, and thus does not lead to held amplification much above 
typical pulsar values). However, the expected variations in the parameter values that 
appear in equation (p|) could cause t sc 1 to range between a few months and a few years. It is 
instructive to compare these values with Atoyan’s (1999) estimate (obtained from a fit to 
the Crab radio data) of the initial spin-down time of the Crab pulsar, £ sd < 30 yr (a factor 
^ 20 smaller than previous estimates that assumed a fixed functional dependence of the 
spin-down torque on fi). The similarity of these estimates is consistent with the possibility 
that the same modeling framework may apply to both plerionic SNRs and SMNS-driven 
bubbles. 


Atoyan (1999) suggested that the initial rotation energy of the Crab pulsar was 
comparable to that of the supernova explosion that produced it, and noted that his inferred 
value of E rot (^ 10 51 ergs) was consistent with independent arguments flChcvalicr 19771) 
that the Crab nebula had originated in a normal Type II supernova event. In the case of 
an SMNS it is, however, unlikely that the explosion energy is as large as the initial rotation 
energy (~ 10 53 ergs), but since the energy (A E Iot ) deposited in the PWB is evidently 
of the same order as E Iot , the supernova ejecta (subscript ej) would be accelerated by 
the expanding bubble and one could obtain an approximate equality between E rot and 
E e j = 0.5M e jU 2 j. For typical ejecta mass ^ 10 M 0 , this would imply n e j ~ 0.1c at t — t sd 
(about an order of magnitude higher than in a typical SNR). This estimate of v e j (which 
agrees with that of VS) is supported by measurements of X-ray emission (e.g., |Piro ct ah| 
2000| ) and absorption (e.g., [Lazzati et al. 20011) features in some GRB sources (see § 13 . 2|) . 


In the supranova scenario, the GRB is associated with the collapse of the SMNS, which 
occurs at a time t sd after the supernova explosion. Unless the explosion takes place within 


14 In view of the large estimated value of ~/ w , we henceforth set (3 W equal to 1 in all numerical expressions. 

15 The total rotational energy of the SMNS is given by E rot = jGM 2 P*/2c, where the parameter j measures 
the stellar angular momentum in units of GMf/c and has values in the range 0.57 — 0.78 for realistic equations 
of state (e.g., Cook et al. 1994; Salgado et al. 1994). 
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a dense molecular cloud, the mass of the swept-up ambient medium will remain negligible 
in comparison with M e j over this time scale and will not affect the SNR dynamics. To 
simplify the treatment, we assume that this is the case.Q The expanding PWB is expected 
to compress the ejecta into a thin shell and accelerate it (e.g., [Reynolds fe Chevalier 198l|) . 
To within factors of order 1, the outer radius of the bubble at time t s d can be approximated 
by the product of v b = v e j(t sd ) times the SMNS spin-down time: 

Rb = v b t sd = 9.5 x 10 m (3 b _iT sd cm , (5) 


where we set v b /c = f3 b — 0.and t sd = r sd yr. To the extent that v b oc (AE rot / M^) 1 / 2 
has nearly the same value in all sources, the magnitude of R b is determined by that of 
t sd . In a similar vein, if the energy AE ro t = 10 53 Ai?53 ergs lost during the SMNS lifetime 
is approximately constant from source to source (AE 53 ~ 1), then t sd can also be used to 
parameterize the SMNS wind power: L w = AE rot /t sd = 3.2 x 10 45 Ai?53/r s d ergs s^ 1 . 


In their original proposal, VS focused on the expected effect of the supranova ejecta and 
SMNS energy release on the baryon content of the environment in which the GRB occurs. 
This was motivated by the general requirement (see, e.g., Piran 1999) that the burst energy 
be channeled into a region with a relatively low number of massive particles in order for the 
outflow to attain the high 10 2 ) Lorentz factors inferred in GRBs. However, this property 
of the GRB outflow is probably determined primarily by the generic properties of the 
central object (e.g., Mcszaros fc Rees 1997 ; Kluzniak fc Ruder man 19UH ; Vlahakis fc KonigI 
2001|) rather than by the matter-sweeping action of the ejecta and SMNS wind. Instead 
of this aspect of the supranova scenario, we emphasize here the favorable consequences of 
the expected delay between the supranova explosion and the GRB event to the creation of 
PWBs in which afterglows with high inferred values of e# and e e could naturally ariseQ 


16 As a further simplification, we neglect the possible incorporation of mass from the ejecta shell into 
the bubble interior through evaporation by the “hot” sliocked-wind material. This effect, which has been 


considered in the study of interstellar bubbles (e.g., Weaver et al. 1977), would be strongly suppressed if the 
magnetic field were strictly transverse to the flow direction, as is assumed in our model. However, even a 
small mean radial field component might lead to a thermal conductivity that is high enough to significantly 
affect the mass budget inside the bubble. 


17 Another potential implication of this delay, which we consider in §k2, involves the interpretation of the 
X-ray spectral features detected in some GRB sources. 
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2.2. Wind-Bubble Structure 


We follow previous treatments of PWB structure ( |Rees fe Gunn 1974| ; [Kennel <£| 
Coroniti 1984| ; [Kmmering fe Chevalier 1987]) in our assumptions about the basic morphology 
of the bubble: we take it to be spherical, with an outer radius Rb, and assume that the 
pulsar wind propagates freely (with a w = const) until it is shocked at a radius R s . Our 
model differs, however, from previous treatments in that we take account of nonthermal 
radiation losses (which could be important during the early phase of the nebula) and we 
do not assume that ideal MHD is necessarily applicable throughout the shocked-wind 
bubble. As has been demonstrated in the previously cited papers, a PWB that expands 
adiabatically with a nonrelativistic speed and that contains a large-scale toroidal magnetic 
held frozen into the matter corresponds to a w ~ /3b 1. Such a model thus cannot 

describe a bubble with j3b -Cl and a w <; 1.Q If a w ~ 1, then the postshock how is 
magnetically dominated from the start. But even if a w -C 1 and the postshock value of 
the fluid-frame magnetic-to-particle pressure ratio ps/p is < 1, this ratio will grow with 
radius r in the bubble and, if radiative cooling is even moderately important and hux 
freezing continues to hold, its value will at some point increase above 1 and could eventually 
become 1. However, as was already argued by Rees & Gunn (1974), a situation in 
which pb significantly exceeds p is unlikely to persist in a real PWB. We therefore adopt 
an alternative formulation and drop the assumption of ideal MHD in the shocked gas at 
the point where the electromagnetic pressure hrst rises above the particle pressure. We 
assume, instead, that beyond that point the electromagnetic pressure in the bubble remains 
in approximate equipartition with the particle pressure]^] For definiteness, we assume that 
the flow obeys ideal MHD within the wind shock and we fix the electromagnetic-to-particle 
pressure ratio 5 = ( S' 2 + B' 2 )/8np in the bubble by setting 


5 = const = max{(5p S , 1} , max{R s , R eq } < r < R b . (6) 

Here £' and B' are, respectively, the fluid-frame electric and magnetic fields, the subscript 
ps denotes postshock quantities, and R eq is the radius where Pb/p first increases to 1 
(assuming ideal-MHD evolution) if h ps < 1. According to this prescription, when the 


18 The presence of a large-scale toroidal magnetic field also implies that the bubble will become elongated 
and will assume a cylindrical, rather than a spherical, shape (see Begelman fe Li 1992|) . We return to the 
question of the bubble morphology in 


19 An equipartition assumption (between the thermal and magnetic pressure components) was previously 
incorporated as a limiting case in the (Newtonian) plerion evolutionary model of Reynolds & Chevalier (1984). 
Possible physical mechanisms for the breakdown of ideal MHD when pb increases to p were discussed by 
Kennel, Gedalin, & Lominadze (1988), Begelman (1998), and Salvati et al. (1998). 
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postshock value of 6 is less than 1 (5 ps < 1 ), the bubble flow starts out maintaining flux 
freezing, but if ps/p increases to 1 before the outer boundary is reached, it switches to a 
nonideal evolution (with S fixed at 1) beyond R cq . If, however, a w is large enough that 
5 ps > 1, then the entire PWB volume is subject to the euipartition constraint, with 5 fixed 
at 5 ps . For a strong shock, the postshock value of 5 can be expressed as a function of the 
wind magnetization parameter a w and speed ft w ^\ Specifically, 


B'J 4 7 ^ s + 1 _ 

8 vTPps ^w^sftps - 1 + (ipsftps/^wftw ) 2 ’ 


(7) 


where 7 ps = (1 — /3 ps ) 1 / 2 is the Lorentz factor of the postshock flow and /3 ps is given by the 
solution of the equation 


ft 


3 

ps 


Au(4 + 5 a w + q2 , A» 


3(& w + 1) 


ftps + 


+ 


ftwOv 


3 (& w + 1) 


= 0 


( 8 ) 


In the limit 7 ^ 1, both <5 ps and f3 ps become functions of a w alone (see GK): equation ([/]) 

simplihes to 


5p S = 


B VS 


^Ipsftps + 1 


8vrp ps 2a- 1 7 2 s /3 2 s - 1 


(9) 


whereas the equation for j3 ps is reduced from a cubic to a quadratic, with the solution 


1 + 2 cr w + [16cr„,(l + a w ) + 1] 1//2 

6(1 + a w ) 


( 10 ) 


For cr^ = 10~ 3 , 10 2 , 0.1, and 1, these expressions yield {/9 ps , <5 ps } = 

{0.33, 0.006}, {0.35, 0.059}, {0.43, 0.53} and {0.73, 4.37}, respectively. Note that, 
for a w < 0.2 (corresponding to (3 ps < 0.5), <5 ps is < 1, so that (by eq. [[]]) 5 = 1. 


We simplify the treatment of the bubble interior by assuming that the flow is purely 
radial, that [in spherical coordinates (r, 6, (ft)] the magnetic held continues to possess only 
a (ft component (B = B 0; see [Bcgclman 1998|) and the electric held only a 9 component 
(£ = £ 6] this follows from the previous two assumptions when ideal MHD is applicable, but 
needs to be postulated when it is not), and that the only nonzero spatial derivatives are in 
the radial direction. Under these assumptions, the particle number, energy, and momentum 


20 In typical applications, the speed of the wind shock is much lower than that of the wind, and even than 
Vb (see Emmering & Chevalier 1987). For the sake of simplicity, we therefore set it equal to zero and identify 


the rest frame of the wind shock with that of the source. 
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conservation equations in the PWB take the form 


+ ;r|:(r 2 7/3n) =0, 


4tt / 


_ 7 A 


;s + 


£B\ 


_ , + AA 

An J ' r 2 dr 


( 


H 2 l3 2 w + 


1 9p _ _t/3A 
r dr c 


(11) 

( 12 ) 

(13) 


where n is the particle number density, w = pc 2 + e + p is the enthalpy density (with e 
being the internal energy density and p the rest-mass density), A is the enrissivity (which, 
like the preceding quantities, is measured in the fluid rest frame), [3 is the radial speed in 
units of c, and 7 is the Lorentz factor. The electric and magnetic held amplitudes in these 
equations are measured in the central-source frame; they are related to their fluid-frame 
counterparts through the Lorentz transformations 


£' = 7 (£ - (3B) , B' = 7 {B - (3£) . 


(14) 


The evolution of £ is governed by Faraday’s law, 


IdB 
c dt 


ld_ 

r dr 


(r£) = 0 . 


(15) 


Given that 7 ^ ^ 1, the shocked gas should be well described by a relativistic equation 
of state 



(16) 


The acoustic speed u ac would be correspondingly high. For example, in the ideal-MHD 
regime, where it is given by [(1/3 + <5/2)/(l + 5/2 )] 1//2 c (representing the phase speed of 
the fast-nragnetosonic wave; e.g., |Konigl 1980|) , v ac /c ~ 0.75 for 5 ~ 1, which we take to be 
large enough in comparison with /7 (^ 0 . 1 ) to justify a steady-state approximation within 
the bubble. We therefore set d/dt = 0 in equations (|TT|)-(|T5|). Equation (|TT| ) then yields 


r r y(3n, -R s 7ps/3ps^ps — G • 


(17) 


The constant C can be evaluated from the shock jump condition 


'7ps/3ps'^ps 7 wfiwH'W - (1®) 

where we have assumed that there is no pair production at the wind shock. Using also 
equation ( 0 ), one obtains 

C _ 

4vr(l + a w ) , y w m e c 3 ' 


( 19 ) 
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Under the steady-state assumption, equation (|15|) implies that 8{r) oc 1/r inside the PWB. 
Normalizing to the value immediately behind the shock, we infer 


£{r) 


£ps-^?ps 


/ Vw PwL w \ 1/ 2 1 
V 1 + &W c ) r 


( 20 ) 


Elimination of the radiative cooling term from equations (|12D and 

(rB) = 0, 


leads to 


dp 2 n dp B - (58 d 
— + 7 fdw— + 
dr dr 


An r dr 

whereas subtraction of y/3 times equation (O) from 7 times equation 


yields 


Id 2 dp 7 (8 — (3B ) d A 

r * l(3w ) _ l/3 JL + n m ; —(rB) = - 

r z dr dr Anr dr c 


( 21 ) 


( 22 ) 


(the entropy equation), where in both cases we took account of the constancy of the 
product r8 inside the bubble. If a w is not -C 1, then most of the bubble volume will be in 
the equipartition regime (h ~ 1, or, equivalently, e# ~ e e ), in which case A will typically 
be dominated by synchrotron radiation.^ To simplify the treatment, we take synchrotron 
emission to be the main radiative cooling process even for low values of a w . Furthermore, 
we assume that at any given location within the bubble the e ± pairs have a monoenergetic 
energy distribution characterized by a random (or “thermal”) Lorentz factor 7 e . The 
latter approximation is appropriate if the postshock gas undergoes significant radiative 
cooling (e.g., |Granot, Piran, fe Sari 2000|) , which, as we discuss in q | 2.3| . may be the case in 
SMNS-driven bubbles. The synchrotron emissivity can then be written in the form 


4 B' 2 

A = -(7 T cn7 e 2 — , (23) 

where <Jt is the Thomson cross section. In view of equation (|l6|), it is then also possible to 
write the particle pressure as 

P = - ^ e nm e c 2 . (24) 

Combining this expression with equation ([lj]) gives 

7 e = Dr 2 ~if3p , (25) 


21 In the equipartition region, synchrotron self-Compton emission cannot exceed the synchrotron radiation 
under any circumstances: it is comparable to the synchrotron emission if the bubble is highly radiative, but 
it remains much smaller if the radiative cooling time is longer than the bubble expansion time. 
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where D 
equation 


3/m e c 2 C = 3 / 7 ps /lp S n ps R 2 m e c 2 . Using equations (p~|), ([T9D , (|24]), and (|25| ) in 
, one can express the radiative cooling term in equation (j 22 |) in the form 


A = G^ 3 p(rB - 13r£) 2 p 2 , 
c 


(26) 


where G = (cr T /2nm e c 2 )D = 6 cr T (l + a w )j w /m e cL w . 


We now consider the term rB that appears in equations (pl|) , (22), and (26). Its form 
depends on whether the flow is in the ideal-MHD or the equipartition regime of the bubble 


interior. The ideal-MHD case corresponds to setting £' — 0 in equation (|14|), which implies 
rB = rS/(5 oc 1/(3 (by eq. [^0|). It is then straigthtforward to obtain from equations 
(pl|) and (^2j) (after also substituting w = Ap from eq. [[T(J) the following pair of coupled, 
first-order, ordinary differential equations for the variables (3 and p, which give the structure 
of the ideal-MHD sector of the PWB: 


d(3 

dr 


Gr 2 £ 2 p 2 + 8 7 2 /3 : 


.2/32 P 


r J 


4 7 4 /?(3/3 2 — l)p — 


3 £ 2 
An (3 


-1 


Rs < r < min{R cq , R b } , (27) 


dp 

dr 


A^ 2 (3 2 p - 


£ 2 


An^ 2 (3 2 


Gr 2 £ 2 p 2 + 87 2 (3 


2o2 P 


n -1 


47T 


£ 2 + A^(3 2 {l-3(3 2 )p 


R,<r< min{R„, R b }, 


(28) 


with £ given by equation 


For the nonideal (equipartition) regime, we combine equations (] 6 j) and (|14|) to obtain a 
quadratic equation for rB , whose relevant root is 


rB — _—_ 

r£l q 1 + pi 


r£ + 


[ 87 t< 57 2 (1 + (3 2 )r 2 p — r 2 £ 2 ] 1/2 
7 2 (1 + (3 2 ) 


(29) 


Using this relation as well as equations (|l 6 l) , (|20|), and (^ 6 |) in equations (^lj) and (|22|) , it is 
once again possible to extract explicit differential equations for (3 and p\ 


d/3_ 

dr 

+ 


1 + 


S(rB eq - (3r£) 


7 2 [(1 + (3 2 )rB eq - 2 (3r£}_ 


8(3 + 


2 5r£ 


^[(l + (3‘ 2 )rB eq -2(3r£] 


p 

r 


GY(3(rB eq - (3r£) 2 p 2 + 
14 7 2 (3 (3 2 - 1) - 5 p+ 


r 2 £B eq — 8n5^ 2 (3r 2 p) [A(3rB eq — (1 + 3(3 2 )r£] 
2n{l + (3 2 ) 7 2 [(1 + (3 2 )rB cq — 2 (3r£]r 2 

max{R s , R eq } < r < R b , 


( 30 ) 
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dp 

dr 

x < 


L , S(rB eq - (5r£) \ - 1 

1 7 2 [(l + /5 2 )r5 eq -2^]J 

(/?r£ — rB eq )(r 2 £B eq — 8716 ^ /3r 2 p) 
27r(l + /? 2 )7 2 [(1 + /3 2 )rB eq — 2/3r£]r 2 
max{R s , R eq } < r < R b , 


47 2 /3p 


(d0\ 25(f3r£ - rB eq ) pi 

\ dr J 7 2 [(1 + P 2 )rB eq — 2(3r£] r j 

(31) 


where the term ( df3/dr ) in equation (|3T| ) is given by the expression (|30|). Equations 
m-m are integrated over their respective validity domains subject to the boundary 
conditions 

P(Ra) = Pps , p(Rs) = Pp S , (32) 

where /3 ps is given by equation (|10|) and the postshock pressure 


Pps 


c r w r YlPln w (R s )7n e c 2 

25 P s7 2 s/3 p 2 s 


(33) 


(with n w (r ) given by eq. [Q) is similarly obtained from the wind-shock jump conditions 
(see GK). The value of R s , where the boundary conditions ([J2]) are imposed, is not known 
a priori and must be determined from an additional constraint. This can be provided by 
requiring global particle conservation: for a bubble considered at time t after the supranova 
explosion, the total number of particles within the radius Rb(t) [which consists of the 
unshocked wind at r < R s (t) and the shocked wind at r > R s (t)] is equal to the total 
number of particles injected by the central neutron star over the time t. The pair injection 
rate at the source is given by 


N 


Lyj 

(1 + cr w )^ w m e C 2 ’ 


(34) 


and hence the total number of particles within Rb at time t is N(t) = Nt. We approximate 
t ~ R b /(3bC , which should be accurate to within a factor of order 1 (for example, 
t = l.bRb/fibC in the case of an adiabatic bubble, with the numerical coefficient decreasing 
in the presence of cooling; see [Reynolds fe Chevalier 19841 ) . The number of particles within 
the volume occupied by the unshocked wind is thus 


JD 

N(r < R S ) = N ' 


p w c ct R b 


whereas the total number of particles within the shocked-wind region is 

r R b „ /3 b N r R b dr 


rUb M, /V r 

N(R S < r < R b ) = / Anr 2 pndr = —— 

Jr 3 R b jb 


R a P 


(35) 


(36) 


where we used equations ®. o. (p~9|) , and (|34|). The solutions obtained in this manner 
will not, in general, be entirely self-consistent, since the bubble structure evolves with time 
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whereas we have assumed a steady state. The wind gas cannot arrive at r = R b at the 
speed Vb after traveling from the origin for the same duration (viz., the age of the bubble) 
as the ejecta that is currently at Rb, given that the wind speed is > Vb everywhere within 
this region and that the ejecta speed was < Vb before it reached Rb . This argument implies 
that, if the particle-conservation condition is imposed, then (3(Rb) will be lower than (rather 
than exactly match) /3b■ This is not a serious inconsistency, since the flow near Rb will 
generally be highly subsonic and therefore can readily adjust to match the speed (u&) of the 
outer boundary. Nevertheless, in order to assess the sensitivity of the results to the choice 
of the imposed constraint, we also solve the system of equations subject to the alternative 
condition (3(Rb ) = /3&. 


2.3. Illustrative Solutions 


The governing equations can be rendered dimensionless by introducing r = r/Rb (with 
r s and r eq denoting the dimensionless counterparts of R s and R cq , respectively) as well as 
B = B/S(Rb ) and p = p/p re f, where 


£ 2 (Rb) _ f 2v w \ ( /3 W L W \ _ f 2a w \ 
Prei 8 tt VI + cr w ) \1 QtxRIc) \1 + (t w ) Pi 


(37) 


(see eq. [^0 []). 22 They then take the form 


d(3 

dr 


+ 4 . 

87 r r 


2 7 4 /3(3/? 2 



r s <r < min{r 5 , 1}, (38) 


dp 

dr 

x 


4 YP 2 p 


rj2p2 r 2 


|^P 2 + 47 2 /3 2 - 
87 t r. 


x 


3 + 2 7 4 /? 2 (l - 3 (3 2 )P 


L r* 


v 13 — 

1 1 J eq 


, r s < r < min{r, 5 , 1} , 

2 y 2 /? + [<5 7 2 (1 + (3 2 )f 2 p — l ] 1 / 2 

7 2 (1 + f3 2 ) 


d(3 

dr 


1 + 


5(rB eq - (3) 


7 2 [(1 + /^ 2 )7 B eq — 2(3} 


3 a 
167T 


- /?) v+ 


(39) 

(40) 


22 The parameter pi is defined so as to isolate the <Tu,-independent part of p le i- it is equal to p le f evaluated 
at a w = 1. 
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+ 


2 fi + 


27 4 [(1 + (3 2 )fB eq — 2/3] 


P 

r 


{[ 7 2 (3/? 2 -l)-5 


p+ 


(f^cg - 8^ 2 fif 2 p)(AfifB eq - 1 - 3/3 2 ) 1 
(1 + /3 2 ) 7 2 [(1 + fi 2 )rB eq — 2/3]r 2 J 
max{r s , r^} < r < 1, 


(41) 


dp 

df 

x < 4 


A , . 

1 7 2 [(l + /3 2 )75 eq -2/?]J 

{fi - fB cq )(fB eq - 8^ 2 fir 2 p) 

(1 + / 3 2 ) 7 2 [(1 + fi 2 )fB eq - 2fi]f 2 

max{r s , ry} < ? < 1. 


7 2 /^P 


/ <#A 2 8(fi-fB eq ) p} 

\ dr ) 7 2 [(1 + /3 2 )f B eq — 2/3] 7 J ’ 

(42) 


(corresponding to eqs. 0-H, respectively), where 


a = 


cttL v 


(1 + cr w )m e c 3 R b 


(&wfiw) 7 w — 


2a 2 


1 + a u 


CL\ . 


(43) 


The parameter cp in equation ( |43| ) was introduced so as to isolate the cr w - independent part 
of the parameter a: it is equal to a evaluated at a w = l.Q As we explicitly demonstrate 
below, ai measures the relative importance of radiative cooling within the bubble. 
Numerically, 1/oq is of the order of the nominal radiative cooling time of the bubble in 
units of Rb/c, and hence the larger the value of ai, the stronger the role that radiative 
cooling plays in determining the bubble structure. In the supranova model, if fi b and A E rot 
are approximately constant from source to source, then oq scales with the bubble age t s ^ at 
the time of the GRB as roughly (see eqs. [Q and 0). 


The boundary conditions given by equation ( |H2| ) are applied at r Sl which, in turn, is 
given either by the relation 


f3 b f 1 dr 


1 - fi b r s Jf s fi(r) 

(see eqs. |54|]-||3(|) or by imposing the condition 

P(f = i) = ft 


= i 


(44) 


(45) 


We solve this system of equations for given choices of the parameter a w and by iterating 
on the value of r s until both the boundary conditions (|32| ) and the constraint ([l4|) or (|45| ) 
are satisfied. 


23 Note that the term in square parentheses in the expression for a can be interpreted as a “bubble 
compactness parameter,” using the terminology often employed in studies of compact astrophysical pair 
configurations. 
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Once the values of /3(r) and p(r) are known, they may be used to obtain the other 
physical quantities of interest. In particular, 


n(r) 

7 e(f) 
'Yw 


_4pi_ 

(1 + a w )(3 w ^ w m e c 2 


1 



r 2 7 (r)/?(r)p(r) . 


(46) 

(47) 


The postshock value of n is given from equations © and © and can, in turn, be used 
with equations ([ 


and ( |33| ) to yield the postshock value of q e , 


7' 


e,ps 


3/3™ {27ps/3ps - [1 - (7ps/3ps/7™/3™) 2 ]} 


lu 


2( 4 7ps/3p S + !)7 P sA 


(48) 


ps 


Furthermore, B(f) is given by l/f(3{r) in the ideal-MHD zone and by equation (|4(]) in the 
equipartition region. We will also find it useful to consider the variable 


^ = 


(S' + B’) 2 

8np 


(49) 


which is equal to Pb/p in the ideal-MHD regime. Note that the ideal-MHD sector of the 
PWB corresponds to the region where 0(f) is < 1 and that, within this sector, -0 = 5. 

Figure |1] shows f s = R s /Rb, the ratio of the wind-shock radius to the outer bubble 
radius, as a function of the cooling parameter a± for the two alternative constraints (N 
conservation and (3(Rb) = (3b) discussed in q[ 2~^ . Results are shown for a w = 1, 0.1, 10~ 2 , 
and 1CT 3 . According to our adopted scalings, Rb oc r S( j oc a[ 1//2 , so Rb decreases with 
increasing Gq. It is, however, seen that when cooling becomes important, the relative 
width of the bubble, 1 — r s , also decreases with increasing a\. This can be attributed 
to the decrease in the internal thermal pressure brought about by the cooling: a lower 
pressure, in turn, requires a shorter length scale to achieve the needed pressure gradient for 
decelerating the flow. For given values of ai and cr w , r s is larger when the terminal speed 
is fixed than when particle conservation is enforced. This can be attributed to the fact 
that, in the former case, the gas speed between r = R s and r = Rb decreases from /3 ps (eq. 
[f0|) to (3b , whereas in the latter case it decreases from f3 ps to (3{Rb ) < (3b (see discussion 
at the end of § |2 . 2| ) : the larger velocity decrement evidently requires a longer deceleration 
length (R b — R s ). Figure |l] also depicts the 6q dependence of r eq = R cq /Rb, the normalized 
equipartition radius (where Pb/p first increases to 1 if <5 ps < 1). It shows that, whereas the 
cr w = 1 solution obeys nonideal MHD throughout the shocked-wind bubble (since <5 ps > 1 in 
this case), the nonideal regime becomes progressively smaller with decreasing a w . This is 
a direct consequence of equation (|9]), which indicates that <5 ps scales approximately linearly 
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with a w : the lower the value of <5 ps , the longer it will take for ps/p in the postshock flow 
to rise above 1 (our condition for the termination of the ideal-MHD regime). The relative 
extent of the ideal-MHD region is larger for the fixed-(3(Rb) solutions, so much so that the 
solutions of this type with the two lowest values of cr w contain no nonideal-MHD zone. This 
can be understood from the systematically higher values of / 3(f ) (see Figs. and [|) and 
correspondingly lower values of B(f) oc 1 /f(3(f) and hence of Pn('f) jp(f) in the fixed-/? (Rb) 
solutions in comparison with their IV-conserving counterparts. 


Figures ^ and |3] exhibit the r dependence of various quantities of interest in PWB 
solutions obtained by imposing the IV-conservation and fixed ~(3(Rb) constraints, respectively, 
for the same 4 values of the wind magnetization parameter as in Figure [T] (spanning 
3 orders of magnitude in a w ) and for 5 values of the cooling parameter (spanning 4 
orders of magnitude in ai). If (3b,-i and A E 55 are both set equal to 1, then the chosen 
values of a\ (= 0.45, 5.04, 45.3, 504, and 4539) correspond to SMNS spin-down times 
t s d = 100, 30, 10, 3 and 1 yr, respectively. Bubbles with Gq near the lower end of this range 
resemble adiabatic PWBs, whereas configurations with ci\ near the upper range are highly 
radiative. Radio pulsars have inferred surface magnetic fields in the range 10 12 - 10 13 G, 
so, by equation (ff[), a variation of roughly two orders of magnitude in the value of t sd is 
naturally expected. 


The upper panels in Figures |2] and |3] show the behavior of (3(f). The dashed lines 
in these panels indicate the postshock (/3 ps ) and outer-boundary (/3b) speeds, whereas the 
dash-dotted curves depict the purely adiabatic (ai = 0) solutions. The displayed results 
confirm that /3(Rb ) < (3b in the IV-conserving solutions. The greatest discrepancy between 
/ 3(Rb ) and the actual speed ((3b) of the outer boundary occurs for a w = 1, in which case 
(3(Rb)/(3b decreases from ~ 0.13 to ~ 0.0063 as a\ increases from 0.1 to 10 3 . We consider 
this to be a tolerable discrepancy, given (as we already noted in § |2.2| ) that the flow near 
r = 1 is highly subsonic. The (3(f) curves further demonstrate that the values of (3 in the 
IV-conserving solutions are lower than those in the corresponding (3(Rb) = (3b solutions also 
for all other values of f where the respective solutions overlap. The lower values of (3 lead 
(on account of the particle flux-conservation relation []lTJ|) to systematically higher values of 
n in the IV-conserving solutions (see the second row of panels in Figs. ^]and |3]). The radial 
profiles of n are nearly flat for low values of the cooling parameter, but when radiative 
effects are important and contribute to the compression, n(f) rises monotonically between 
f s and f — 1 with a slope that is steeper (particularly in the nonideal-MHD regime) the 
larger the value of a\. 


The third row of panels in Figures |2] and |3] displays r ) e (f)/lw The dashed line marks 
the postshock value of this quantity, and the dash-dotted curves again represent the purely 
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adiabatic case. The random Lorentz factor declines monotonically with r in the a w = 1 and 
0.1 solutions, but at lower values of a w (for which the electromagnetic pressure contribution 
is small) and so long as cooling remains relatively unimportant (ai 10 ), 7 e (r) initially 
increases behind the wind shock to make it possible to attain the necessary total pressure. 
The minimum value of 7 e (r) (reached at f = 1) decreases (as expected) with increasing a± 
and is lower for the TV-conserving solutions than for the &xed-/3(Rb) ones. This difference 
can be understood from the fact that the particle density near the outer boundary is lower 
in the latter case [corresponding to a higher value of /3(f = 1)], so y e must be larger to bring 
the pressure up to its requisite value. 

Of particular relevance to the evolution of GRB afterglows is the behavior of the 
thermal and electromagnetic pressures. The radial profiles of p (normalized by pi) are shown 
in the fourth row of panels in Figures and [3]. It is seen that the behavior of p(f) varies 
with the parameter choices and that the details depend on the nature of the constraint 
imposed on the solution. For a w = 1, the curves decline monotonically when cooling is 
relatively unimportant but increase monotonically at high values of a\. For a w = 0.1, the 
curves decrease monotonically for all plotted values of the cooling parametr, whereas for 
a w = 10 -3 in the fixed-/3(/4) solution they increase monotonically for all exhibited values 
of a\. In the remaining cases, the curves increase with r near the inner boundary of the 
bubble and decrease near its outer boundary. Since p oc 7 e n, one can understand the shape 
of the curves by comparing them with the corresponding curves in the second and third 
rows of panels. One finds that, quite generally, the behavior of p is dominated by that of 7 e 
for low values of ai, but that the influence of the density variations becomes progressively 
more important as the cooling parameter increases. 

The effect of the electromagnetic fields can be inferred from the behavior of the variable 
0 and of the ratio S' / B' , shown in the bottom two rows of panels in Figures ^ and |3j. 
Note that one can write 0 = 5 + S'B'/Anp (see eq. [[49]). Immediately behind the wind 
shock (in which, by construction, S' = 0), 0 is equal to 5 ps (given by eq. 0). If 5 ps > 1 
then the nonideal-MHD regime starts right there: beyond that point, 6 remains fixed at its 
postshock value but 0 and S' / B' increase with r. If 5 ps < 1 , then the evolution proceeds 
under ideal-MHD conditions (S' = 0), with 0 continuing to coincide with 5 (which in 
this regime equals Pb/p) up to the point where it reaches 1. This point marks the end 
of the ideal-MHD regime and corresponds to f eq . Beyond r eq , 0 continues to increase 
monotonically with f (and now so does also S'/B'), but 5 remains fixed at 1 (see eq. ||). 
The figures show that, as expected, the electromagnetic contribution to the total pressure 
becomes progressively larger as a w increases, and they further demonstrate that S'/B' 
exhibits a similar trend. It is also seen that cooling enhances the relative importance of the 
electromagnetic fields, which can be understood from the fact that it reduces the magnitude 
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of the thermal pressure component. The overall behavior of ?/; and £' / B' does not appear to 
depend strongly on the choice of constraint under which the solution is obtained, although 
the corresponding curves in the two figures differ in their details (see also the related 
discussion in connection with Fig. |T] above). As we show in §!0, the variable -0 plays a key 
role in the modeling of relativistic shocks that propagate inside a PWB. 


3. Implications to GRB Afterglows 
3.1. Emission-Region Parameters 

Pulsar wind-inflated bubbles, such as those predicted to arise prior to the onset of 
the high-energy burst in the supranova scenario, provide an optimal environment for GRB 
afterglows since they naturally yield high electron and magnetic energy fractions (e e and €b) 
behind the propagating shock wave that gives rise to the afterglow emission. High values of 
e e are expected from the fact that relativistic pulsar-type winds are likely dominated by an 
electron-positron component, whereas significant values of eg should naturally occur if the 
winds are characterized by a high magnetization parameter. 

The observationally inferred values of e e and €b are derived from spectral fits that are 
based on the standard model assumptions of a “cold,” weakly magnetized proton-electron 
preshock medium. We now consider what would be the “equivalent” values (which we 
denote by the subscript equiv) that one would derive if the afterglow-emitting shock 
propagated instead inside a PWB. The postshock quantities can be determined from the 
appropriate generalizations of the expressions presented in §2.2| . taking account of the fact 
that the preshock gas is now “hot” (with w = 4p; eq. 0) rather than “cold” and that it 
may contain a nonzero comoving electric field. One imposes the continuity of the energy 
flux y 2 (3w + £B / 47 T, momentum flux 7 2 /3 2 w + p+ (£ 2 + B 2 )/8n, and tangential electric field 
£ in the frame of the shock. The possible presence of a nonzero £' requires the specification 
of an additional shock jump condition, which we take to be the conservation of magnetic 
flux during the fluid’s transit through the shock. On the assumption that the field is 
transverse to the shock propagation direction, this condition translates into the requirement 
that /3B be continuous in the shock frame. Combining these constraints, one obtains a 
cubic equation for the postshock (subscript 2) flow speed (measured in the shock frame) 
that is akin to the expression (j^). In the limit of an ultrarelativistic shock, the latter again 
reduces to a quadratic equation, whose relevant solution is 

_ 1 + i/j + [(1 + iJj) 2 + 3(2 + ^)V’] 1/2 
3(2W0) 


( 50 ) 
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where ip is defined by equation (f49[). Equation ( |50| ) reproduces equation (III]) if one 
substitutes 2a w for ip: this is consistent with the fact that ip /2 becomes equal to the 
magnetization parameter a = £ B / Att^ 2 f3w of the preshock medium in the limit f3 —» 1 and 
w —> 4 


By substituting £ = 'y(£' + f3B') and B = 7 (B' + /3£') (see eq. |b|]) into the expressions 
for the energy and momentum fluxes in the shock frame, one finds that, if the shock is 
highly relativistic (so that the upstream (5 is close to 1 in the frame of the shock), the 
upstream fluxes have the same forms as in a purely hydrodynamic shock if w in the latter 
is replaced by w + ( B’ + £’) 2 /47T. Shock models of GRB afterglows traditionally infer an 
ambient gas density by assuming that a hydrodynamic shock propagates into a standard 
ISM or stellar-wind environment with an enthalpy density w = niim p c 2 , where m p is the 
proton mass. This motivates us to define the “equivalent” hydrogen number density 


n H ,equiv 


1 r ( b ' + £') 2 

- 9 W - A - 

m p c z [ 47 T 


(51) 


which under the assumption of a relativistic equation of state can be written as 
^H,equiv = 4(1 + ip/2)p/m p c 2 . This quantity is plotted in the top row of panels in 
Figures |] and [5] as a function of r for each of the model PWBs presented in § [2T3| . 

These figures also show (in the second row of panels) plots of the radial dependence 
of k = — dlognH,equiv/dlogr, the effective power-law index of the equivalent hydrogen 
density distribution. Our model assumption of an abrupt transition between the ideal- and 
nonideal-MHD regimes at r eq introduces an unphysical discontinuity in the value of k at 
this point: the displayed curves have been smoothed at this location by interpolation across 
the discontinuity. 

To simplify the discussion, we restrict attention to the three synchrotron-spectrum 
characteristics considered by Sari et al. (1998, hereafter SPN), namely, the break frequencies 
u m and u c and the peak flux Ph imax ; we refer the reader to that paper for the definition 
of these quantities and for the derivation of the standard expressions for emission by a 


24 An indirect measure of the value of the magnetization parameter just ahead of the afterglow-producing 
shock is possibly provided by the power-law index p of the synchrotron-emitting particle energy distribution, 
which can be deduced from the shape of the observed spectrum (e.g., pari, Piran, fe Narayan 1998| ). Kirk et 
al. (2000) argued that an ultrarelativistic, unmagnetized shock that accelerates particles in a “cold” medium 
through the first-order Fermi process produces a “universal” power law of index p ss 2.2, and that this value 
increases with the preshock magnetization (so that, for example, p ~ 2.3 for er = 0.01, assuming £’ = 0). 
The often-quoted “canonical” value of p for GRB afterglows is 2.5, although in some sources a significantly 
higher value has been inferred (e.g., Huang, Dai, fe Lu 2000 ; Panaitescu fe Kumar 2002 —note, however, 
that the latter reference also lists a few sources in which p < 2.2 has been deduced). 
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spherical shock in which there is no pair production. In the interest of simplicity, we also 
assume that the equivalent hydrogen number density inside the bubble is roughly constant; 
as the U-H.equiv pl°ts i n Figures [3] and [5] demonstrate, this approximation is usually adequate 
over the bulk of the bubble volume, especially when the cooling is not too strong. We 
distinguish between two cases: weakly cooling PWBs (corresponding to cooling parameters 
ai <; 10 2 , or, for our fiducial values, r sd ^ 10), whose radial widths A R b = ( R b — R s ) are of 
the order of R b , and strongly cooling PWBs (cq S> 10 2 , r sd <C 10), for which AR b /R b -C 1. 
In the weakly cooling case, one can approximate the volume of the shocked bubble gas by 
that of the sphere that is bounded by the shock. 

In the standard case of a uniform ambient medium and a slow-cooling (adiabatic) shock, 
one can express the two break frequencies and the peak flux in terms of the shock energy 
E, the ambient density nn, the observed time t, as well as e e , e#, and the distance to the 
source (see eq. [11] in SPN). In particular, u m oc e 2 e]/ 2 F ,1//2 f _3 ' /2 , u c cx , 

and .F^max oc e^ 2 En^ 2 . In the case of a shock propagating inside a weakly cooling PWB, 
it turns out that the above expressions for u m and v c are reproduced if %, e e , and e# are 
everywhere replaced by nn,equiv ( ec f [0); 


L e,equiv 



and 


B. 


/ 2 


e B, 


equiv 


327r7| 1 7r H ,e q uiv?«pC 2 



^e'Je > 

(52) 

f. 

(53) 


respectively, where 72 is the Lorentz factor corresponding to /3 2 , 721 is the Lorentz factor of 
the postshock fluid in the stationary frame (evaluated again on the assumption that the 
afterglow-emitting shock is highly relativistic), and y e is given by equation ([47]). In deriving 
equation (|52|), we have made use of equation ( |24| ) and of the continuity of the particle 
flux 7 (3n in the shock frame. The expression for F„. max is reproduced by making similar 
substitutions and then multiplying by the factor 


F t 


correct 


n 


nn. 


equiv 


3 m p 

4 m e 



(54) 


The functions e e ,eqmv(h), eB !e qmv(h), and ^correct (h) are plotted in the third-through-fifth 
rows, respectively, of Figures £| and [|. with the dashed line in each panel marking the 
value of the respective quantity immediately behind the pulsar-wind shock. Although the 
underlying expressions were obtained under the assumption of a weakly cooling bubble, 
it turns out (see next paragraph) that they continue to apply also in the strongly cooling 
case. The values of e eiequiv were calculated by setting e e = 1 in equation (|52|) — consistent 
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with our PWB model approximation of a pure-e 1 * 1 wind. Inasmuch as F correct does not differ 
from 1 by more than a factor of a few over most of the explored parameter range, it can be 
concluded that the standard expressions will remain approximately applicable if one simply 
replaces tih, e e , and e# by their “equivalent” counterparts. The derived values of F correct 
also verify that riH, e quiv is usually much greater than (■ m e /m p )n , consistent with the PWB 
model assumption of a “hot” equation of state. 

To derive the corresponding expressions for a strongly cooling PWB, one can 
approximate the bubble as a thin shell of radius ~ Rb . The bubble volume traversed 
by a shock that is located at a distance x from the inner radius of the bubble is then 
V (x) ~ 47 tR%x. Relating x to the observed time t and to the Lorentz factor T of the shocked 
gas by x ~ 4T 2 t, following SPN,Q and setting E ~ VT 2 n^ equ i v m p c 2 in the adiabatic-shock 
case, one finds that the standard expressions (eq. [11] in SPN) continue to apply if one 
makes the aforementioned substitutions for e e , eg, and tt-h, and, in addition, multiplies the 
expressions for u m and u c by A correct and l/A correct , respectively, where 


( AEt \ 1/2 

correct — I ~TZ 7ZT I 

y 177TC772p72H ,equiv Rb J 


q cm w}l 2 + x l 2 r) ~ 1 l 2 /?-2 
o.oan/ 52 i d n Hequiv rt b l7 


(55) 


Here F 52 = (F/10 52 ergs), t d is the observed time in units of days, and Rb,n = (Rb/ 10 1 ' cm 
(see eq. [S]]). The expression for the flux correction factor remains the same as in the 
weakly-cooling-PWB case; the displayed plots indicate, however, that F umax undergoes 
large variations across the bubble when the particle-conservation constraint is imposed 
(although not when the terminal speed is fixed). Although the factor A correct alters the 
parameter dependences of the break frequencies (specifically, v m oc Et^ l n u l J^ uiv and 
u c oc £ ,_1 i _1 nH ) 1 e / ^ 1 iv for an adiabatic shock in a strongly cooling bubble), its numerical value 
will not be large for typical afterglow parameters. This seems to suggest that the standard 
expressions should provide adequate estimates of the source parameters in this case, too, 
but we caution that the adopted approximation of an effectively uniform ambient medium 
becomes questionable at large values of aiQ 


25 In a more precise treatment, one obtains t for radiation emitted along the line of sight to the center 
from the differential equation cdt/dr = l/2F 2 h , where T s h is the Lorentz factor of the shock (see Sari 1997 ). 
The solution in this case is t = a;/4r 2 h c + R s /2TqC (where To is the initial Lorentz factor of the outflow and 
Tsh ~ [(1 + /? 2 ) /( 1 — /? 2 )]r for r sh 1), which shows that the approximation used in the text is only valid 

for x »(r/r 0 ) 2 R s . 


26 In the conventional interpretation, the forward shock is expected to be fast-cooling (and, if e e is close 
to 1, also radiative) during the early phase of the afterglow evolution. However, in the case of a PWB 
environment characterized by a w <; 1, the shock can be only partially radiative since a significant fraction of 
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The above considerations suggest that, as a rough check of the compatibility of the 
PWB model with observations, one can examine the consistency of the predicted values of 
e e ,equiv, e,B,equiv, and nn.equiv with the values of e e , e#, and % that are inferred from the 
spectral data by using the standard ISM model. The results shown in Figure £| indicate 
that, if a w is not <C 1, then e e) eq U iv typically lies in the range ~ 0.1 — 1. As we noted 
in §1, such comparatively large values have been inferred for the corresponding standard 
parameter e e even before a large body of data became available, based on emission-efficiency 
considerations as well as on some early model fits. These inferences have been supported 
by more recent analyses of the accumulating data on afterglows, which have even led to 
the suggestion that e e may have a “universal” value ~ 0.3 ([Freedman fe Waxman 200 1| ; see 
also [Huang, Dai, fc Lu 2000| and [Panaitescu fc Kumar 2002|) . The predicted magnitudes of 


-e,equiv 


are larger for smaller values of a w and when the fixed-/3(i4) constraint is imposed 
(see Fig. |5j). However, these estimates could be lowered in real sources by the admixture 
of baryons into the bubble interior — either by injection at the source or by evaporation 
from the bounding supranova-ejecta shell (see footnote 16). Turning next to the magnetic 
energy-density parameter e^equiv, we see from Figures [| and [5] that values in the range <; 1 
to <; ICG 2 are predicted by our model as a w decreases from 1 to 10~ 3 . It is noteworthy 
that this range is consistent with the values inferred in the standard ISM picture for a 
source like GRB 970508 (which, as we discussed in §1, have posed a challenge for the 
conventional scenario) as well as with the mean values of e# inferred for the afterglow 
sample of Panaitescu & Kumar (2001). In this case one can again expect a reduction in 
the estimated parameter value as a result of the admixture of baryons (which would reduce 
the cooling and the associated magnetic field compression in the bubble; see GK), although 
other factors may also contribute to a lowering of e#: for example, the afterglow-emitting 
shock may not be transverse, and the electromagnetic-to-thermal pressure ratio in the 
diffusive regions of the PWB might be lower than the equipartition value adopted in our 
model. 

Finally, our derived values of ?tH,eqmv are compatible with the observationally inferred 
preshock particle densities. The values of rip estimated in the literature under the 
assumption of a uniform, “cold” ambient medium typically span the range ~ 0.1 — 50 cm~ 3 
(e.g., IPanaitescu fe Kumar 20021) , although (as noted in the last-cited reference) there 


are examples of sources where a density < 10 2 cm 3 is implied. There have also been 


the shock kinetic energy is converted into postshock magnetic energy, which is not subject to radiative losses. 
Although one can in principle derive the appropriate expressions also for this situation [as was done, e.g., 
by Bottcher & Dernier (2000) in the standard case], we do not consider fast-cooling shocks here since the 
results for partially radiative shocks are more cumbersome and would unduly complicate the presentation. 
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suggestions in the literature that some afterglows originated in a medium with a density 
> 10 2 cm -3 . As is seen from Figures [4] and [j]. our model can in principle account for 
all of the inferred values: the typical densities are reproduced by PWBs with r s d in 
the range ~ 3 — 30, whereas “outlying” low and high inferred densities correspond to 
lower (respectively, higher) values of the cooling parameter cii. The basic trend is for 
bubbles with more radiative cooling to be characterized by higher values of the equivalent 
density p 7 ] this follows directly from our adopted parameterization (see § [2.1|) , which implies 
^H,equiv oc pi oc r sd oc a/ . It is, furthermore, seen that the effective-density predictions are 
remarkably insensitive to the choice of the wind magnetization parameter a w and of the 
specific constraint imposed on the solution. This robustness can be traced to the fact that 
?ZH,equiv basically measures the energy density in the bubble, which, for a given choice of 
(and thus of Rb oc aj~ ^ 2 ), is essentially determined by the wind ram pressure at a distance 
rs^ Rb from the center. 


Another attractive feature of the PWB scenario is that it naturally gives rise to radial 
profiles of nn.equiv that, depending on the cooling parameter a± and the location within 
the bubble (see the plots of 2 — k in Figs. |4] and |5|), may resemble a uniform medium 
(constant-nn ISM or interstellar cloud) or a stellar wind (up oc r~ 2 ; but note that k strictly 
remains ^ 1 in these solutions). Both types of behavior have, in fact, been inferred in 
afterglow sources (e.g., |Chcvalicr fc hi 2000| ; [frail ct al. 2000| ; |Halpcrn ct al. 20001) . The 
unique aspect of the radial distribution of nn, e quiv in this picture is that it spans a range 
of effective power-law indices k that can vary from source to source, and, moreover, that 
the value of k appropriate to any given afterglow is predicted to change with time as 
the afterglow-emitting shock propagates within the bubble. This leads to a more flexible 
modeling framework for the afterglow evolution and can naturally accommodate cases 
where a value of k that is intermediate between those of a uniform ISM and a stellar wind 
could best fit the observations (see, e.g., [Livio fe Waxman 2000| ). It also explains why 
afterglows associated with star-forming regions need not show evidence of a stellar-wind 
environment (as expected when the GRB progenitor is a massive star; in view of the derived 
magnitudes of n H , e quivj this model also makes it possible to understand how a source with 
such a progenitor could produce an afterglow with an implied value of % that was much 
lower than the typical ambient density near massive stars). In addition, high values of 
^H,equiv in this picture are not subject to the objection (e.g., |Halpcrn ct al. 2000|) that 


27 In particular, afterglows with inferred preshock densities above ~ 10 cm -3 are expected in this picture 
to arise in highly radiative PWBs, suggesting that such sources may be the most promising candidates 
for testing the predicted departures from the standard spectral scaling relations (which, according to our 
preceding arguments, should be most pronounced in rapidly cooling bubbles). 
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they will necessarily give rise to excess extinction (although it is also conceivable that dust 
destruction by the optical-UV and X-ray radiation from the GRB outflow could reduce any 
preexisting extinction toward the source; see |Waxman fc Draine 20(J0| and [bruchter, Krolik, 
|fc Rhoads 200 1|). As is seen from Figures Q and [^, the predicted r2,H, e quiv(' r ) distributions 


exhibit progressively steeper declines as the outer boundary of the bubble is approached. 
This suggests that the later phases of the evolution of any given afterglow would be more 
likely to exhibit signatures of a stellar-wind environment. In all the cr w > 0.1 bubbles, this 
wind-like behavior becomes more pronounced the lower the value of the cooling parameter 
(i\. Since (as noted above) the value of n H , e quiv also exhibits a systematic dependence on 
this parameter (it decreases with decreasing ai), one may expect afterglows with higher 
inferred ambient densities to preferentially indicate a uniform ISM-likc environment if they 
originate in such bubbles. 


The inferred radii of afterglow shocks typically lie between ^ 10 1 ' cm and <; 10 1S cm 


ps 


(e.g., [Piran 1999| ; |Chcvalicr fc Li 2000| ). These values are consistent with the upper limit on 
the bubble’s outer radius (eq. |J) for supranova-GRB time delays of ^ 1 yr to <; 10 yr. We 
can check on whether typical afterglow-emitting shocks will still be relativistic by the time 
they reach the outer edge of the bubble at Ri by solving the adiabatic evolution equation 


[r 2 (r) - 1 ]M eq (r) + [T(r) - 1 ]M 0 = (T 0 - 1)M 0 , 


JL < r < Rh 


(56) 


where 


M eq (r) 


— / 47 tR ^H^equiv V, 

J Rs 


r {R)m p dR (57) 

and M 0 = E/T 0 c 2 (e.g., van Paradijs et al. 2000). The bottom panels in Figures [| and [5] 


show the results for the PWB solutions presented in § |2.3| using a representative value of To 
and two plausible values of E. It is seen that, in all cases, the GRB outflow decelerates 
rapidly after entering the bubble, and in weakly cooling PWBs the Lorentz factor of the 
afterglow-emitting gas is at most a few by the time the shock reaches R^ (and is effectively 
nonrelativistic for the E = 10 52 ergs solutions). Only in the case of an energetic shock and 
a strongly cooling bubble is r(/4) appreciable (but even then it remains ^ 10). It is worth 
bearing in mind, however, that, if the outflow is collimated with a small opening half-angle 
9j, then it will start to strongly decelerate due to lateral spreading when its Lorentz factor 
decreases to ~ 1/9j (e.g., |Sari, Piran, fc Halpcrn 1999| ; [Rhoads 1999|) , so that even the 
more energetic shocks could become nonrelativistic while they are still inside the PWB.^j 
These results are again quite insensitive to the value of a w and to the choice of the imposed 


28 Mass loading of the bubble by evaporation of the ejecta shell will further contribute to the lowering of 


r. 
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constraint; this is not surprising in view of the fact that they depend mostly on the behavior 
°f ^H,equiv(f9, which exhibits a similar trend. 

A GRB shock that reaches the supranova ejecta shell at r = /?/, with a Lorentz factor 
> 1 would be rapidly decelerated to subrelativistic speeds since the rest-mass energy of 
the shell (~ 2 x 10 55 M e j ilo ergs, where M e j jl0 = M e j/1OM 0 ) is in most cases much greater 
than the (equivalent isotropic) shock energy E. The spectral characteristics of the forward 
shock after it enters the shell could be evaluated once the dynamical evolution of the shock 
is calculated.^ Besides the anticipated alterations in the spectral scaling laws, one may 
expect the numerical values of the various physical parameters to undergo dramatic changes 
as the shock moves from the interior of the bubble to the ejecta shell: in particular, n h 
would likely increase by several orders of magnitude, whereas e e and e# would probably 
decrease significantly. In addition, the sudden deceleration would drive a relativistic reverse 
shock into the GRB outflow, whose emission may have an important effect. The overall 
outcome is likely to be a discontinuous change in the shape and evolution of the observed 
spectrum. Ramirez-Ruiz et ah (2001) modeled a somewhat similar situation that may 
arise when a GRB shock that propagates in a stellar wind encounters a density bump. 
They suggested that an encounter of this type could induce a brightening and reddening 
of the afterglow spectrum and might explain observations of such a behavior in several 
sources. The situation considered by Ramirez-Ruiz et al. (2001) differs, however, from a 
PWB-SNR transition in that the density contrast as well as the jumps in e e and e# (which 
were assumed to be negligible in the “bump in a wind” model) would typically be much 
larger in the latter case. It would thus be interesting to carry out a detailed investigation 
of the observational implications of the shock encounter with a dense shell in the explicit 
context of the PWB model. 


3.2. Interpretation of X-Ray Features in the Supranova/PWB Model 


The SNR shell bounding the PWB could also manifest itself by imprinting X-ray 
features on the GRB afterglow spectrum. Indeed, recent detections of such features in 
several GRB sources have been argued to provide strong support for the supranova scenario 
(e.g., [Lazzati, Gampana, fc Ghiscllini 1999| ; [Piro et al. 2000| ; |V ietri et al. 200l| [Amati 
|ct al. 2000| ; Lazzati et ah 2001 ; Bottcher, Fryer, fe Bernier 2002| ). To date, four GRB 


29 Although the behavior of the shock in both the highly relativistic and the Newtonian limits had been 
considered in the literature, so far there has been no published treatment of the transition between these 
two regimes. 
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sources (GRB 970508, GRB 970828, GRB 991216, GRB 000214), observed -8-40 hr 
after the burst, showed emission features in their postburst X-ray spectrum, and one source 
(GRB 990705) exhibited an absorption feature that disappeared 13 s after the onset of 
the burst. These features most likely represent Fe K a lines or an iron K edge, and their 
detection implies that a large quantity 0.1 M 0 ) of pure iron is located in the vicinity 
(r ^ 10 16 cm) of the GRB source. Such a large iron mass is most naturally produced in a 
supernova explosion, and the inferred distance of the absorber indicates that the supernova 
event preceded the GRB by at least several months, as expected in the supranova picture. 
The association with a supranova is further strengthened by the argument flVietri et al.| 
2001|) that the abundance of 56 Fe (the product of the radioactive decay of 56 Ni and 56 Co) in 
supernova ejecta is not expected to become significant until — 10 2 days after the explosion, 
during which time the ejected gas in a source like GRB 991216 (in which the observed line 
width is consistent with an outflow speed — 0.1 c; [Piro et ah 2000|) would have traveled to 
a distance ^ 10 16 cm from the origin. 

We now proceed to discuss how the observed X-ray features can be interpreted in the 
context of the supranova scenario, and we consider the implications of this interpretation to 
the PWB afterglow model presented in this paper. We concentrate on the specific example 
of GRB 991216, which allows us to capitalize on the analysis already carried out on this 
object by Piro et al. (2000) and Vietri et al. (2001); our interpretation does, however, differ 
in its details from the model favored by the latter authors. We approximate the ejecta as a 
thin spherical shell of radius R c j and density n e j. Although the ejecta of a supernova that is 
not associated with a pulsar may be expected to fill the volume into which it expands, in 
the case of an inflating PWB the ejecta will be swept up and compressed into a dense shell 
(e.g., |Ghcvalicr 1977|) . The acceleration of this shell by the lower-density bubble gas would 


subject it to a Rayleigh-Taylor instability, which could lead to clumping (see, e.g., [Jun 
|1998|) . As we argue below, such clumping is consistent with the data for GRB 991216.Q 

We assume that the emission is induced by continuum irradiation from the central 
region that commences around the time of the burst but is not necessarily confined to the 
solid angle of the GRB outflow. The part of the shell that is observable to us at time t is 
limited by light-travel effects, so that, for a source observed up to time t max , the solid angle 
Afi from which Fe emission is received is given by 


AfT 1 — cos 9 m ax ct 


max 


47T 


2(1 + z)R e j 


1.1 x 10 15 cm 


R, 


(58) 


j ej 


30 The presence of a strongly clumped shell was already inferred by Lazzati et al. (2000) in GRB 990705 
from their analysis of the X-ray absorption feature in that source. 
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where the angle 9 is measured with respect to the line of sight to the origin, and where we 
substituted numerical values appropriate to GRB 991216 (redshift z = 1.02, £ max = 40.4 hr). 
Piro et al. (2000) identified the X-ray feature in GRB 991216 as an Fe XXVI Ha line 
(rest energy 6.97 keV) with a FWHM (as quoted in [Lazzati ct al. 200~T|) of ~ 0.15 c. 
Since we attribute the emission to material that moves toward the observer with a speed 
of that order, we favor an identification with a lower-energy line, specifically Fe XXV 
Hect (rest energy 6.7 keV), but our results are not sensitive to this choice.^] Based 
on the photoionization models of Kallman & McCray (1982), Fe XXV is the dominant 
ion when the ionization parameter £ = L,/n^ T? 2 (where Lj is the ionizing continuum 
luminosity) lies in the range log£ ~ 2.7 — 3.2. Using L* = AnD 2 F x = L,^ 


10 45 ergs s 1 , 


with F x = 2.3 x 10 12 ergs cm 2 s 1 and D — 4.7 Gpc ( [Piro ct al. 200C| ), we thus infer 


n ej /?4 = 6.1 x 10 42 (L ii45 /6.1)(£/10 3 ) 


3 \ —1 


(59) 


The observed line luminosity corresponds to N FBj 52 = (lV Fe /10 52 photons s 4 ) = 8 ([Piro et 
al. 2000|), and we can write N Fe = ( AQ/4-7T ) M Fe / 5 6m p t rec , where M Fe = 0.1Mp ei o.i M 0 is 


the total iron mass in the shell and t rec ~ 4 x 10 9 Tg- 6 Z _2 rig 1 = 2.8 x 10 1 °T'g- 6 ri“ 1 s is the 
recombination time for a Z = 24 ion (with T e and n e being the electron temperature and 
number density, respectively, and Z the ion charge). The expression for f rec is valid in 
the range T e ~ 10 2 — 10 6 K ( [Lazzati et ah 2001|) , and photoionization models imply that 
T 6 = T e /10 6 K « 1 for log£ sa 3 ([Kallman fc McCray 1982|) . Approximating n e ~ n e j, we 
obtain 

n ej = 1.0 x 10 9 (47r/AO)T 6 a6 (7V Ffe)52 /8)M Ffe ; (U cm" 3 . (60) 


Substituting equation (|6D[) into equation ( |59| ) gives 

Re j = 7.6 x 10 16 (AH/47r)°- 5 M 9 e 5 0 . 1 (iV Fei52 /8)- 0 - 5 (L i ,45/6.1) 0 - 5 T 6 - 0 ' 3 (£/10 3 )- 0 ' 5 cm . (61) 

Combining equations ([H]) and (|58|) then yields 

Rej = 1.9 x 10 16 M^ o . 1 (A Fe , 52 /8)- 1 / 3 (L i ,4 5 /6.1) 1 / 3 T 6 - 1/5 (£/10 3 )- 1 / 3 cm (62) 


and (for the given fiducial values) 6 
relation (!>9|) in turn implies 


~ 28°. Substituting the estimate (^2]) into the 
n ej = 1.8 x 10 10 M Fe 2 Q 3 1 (iV Fe! 5 2 /8) 2/3 (Li i 45/6.1) 1/3 T 6 2/5 (£/10 3 ) _1 ^ 3 cm" 3 . (63) 


31 An independent argument for an identification of the X-ray feature in a source like GRB 991216 with 
the Fe XXV Hea line was presented by Ballantyne & Ramirez-Ruiz (2001), who demonstrated that an Fe 
XXVI Ha line is unlikely to be observed because of the removal of photons from the line core by Compton 
scattering. 
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If the shell expands with a speed u e j ~ 0.1c (see m, then its age when it reaches the 
radius given by equation ( |62| ) is 

i age « 72 Mpg Q 1 (lV Fei 5 2 /8) _1 / 3 (Lj i 45/6.1) 1/3 T 6 ~ 1 ^ 5 (^/10 3 ) _1 / 3 (u e j/0.1c) _1 days . (64) 

This value is consistent with the time required for the bulk of the ejected radioactive 56 Ni 
to decay into 56 Fe. 

The electron column density in the X-ray emitting portion of the shell is given by 
N e = Mej/AirRl-fA.He'nT'p, where /a is the covering factor of the shell and /r e , the electron 
molecular weight, is 2 for hydrogeu-free ejecta. Assuming /a ~ 1 and using the estimate 
(|62l), the Thomson optical depth of the shell is inferred to be 

r T = 0.9M ejil0 Mp e 2 0 /3 1 (iV F e,52/8) 2/3 (L i) 45/6.1)- 2/3 T 6 2/5 (e/10 3 ) 2/3 (/i e /2)- 1 . (65) 


For these fiducial values, the thickness of a homogeneous shell would be ~ 8 x 10 13 cm, 
which is consistently <C R e y It is, however, more likely that this nominal thickness 


corresponds to the size of a clump in a shell with a small volume filling factor (see [Lazzati 
|et al. 200T|) . In fact, a high degree of clumping is also indicated by the requirement that 
the line photons reach the observer without undrgoing excessive Compton broadening in 
the shell. The photoionization optical depth of the iron ions in the shell is similarly inferred 
to be 

r Fe = 4.4 Mpgo. 1 (77 Fei 5 2 /8) 2 ^ 3 (Li ) 4 5 /6.1) _2/3 T 6 2/5 (^/10 3 ) 2 ^ 3 (?7/0.5) , (66) 


where t] is the relative abundance of the Fe XXV ion (e.g., [Kallman fe McCray 1982]) and 
where we used a Fe xv ~ 2.0 x 10“ 20 


cm 


T-g., 


Krolik &; Kallman 1987| ). (Our fiducial 
mass ratio M Fe /M e j = 0.01 corresponds to an iron abundance that is ~ 5.6 times the 
solar value.) The estimated values of tt (<; 1) and r Fe (a few) are optimal for producing 
high-equivalent-width iron lines through reflection (e.g., |Weth et al. 20001 ; |Vietri et al.| 
|2001|) . Since, in this picture, r Fe oc 1 /T? A , the efficiency of producing detectable emission 
lines would typically be low for shells with radii much in excess of ~ 10 16 cm (eq. [|62| ). 


The most natural way of relating the above scenario to the PWB model is to identify 
fage with t s d and R e j with Rb. However, such a straightforward identification is problematic 
in that the magnitude of R e j that is inferred from the X-ray emission-line observations 
(<; 10 16 cm; eq. |6l]]) is at least an order of magnitude smaller than the lower limit on 
Rb typically implied by the afterglow data. In particular, in the case of GRB 991216, 
the optical light curve showed evidence for steepening (which was attributed to shock 
deceleration triggered by the lateral spreading of a jet) starting about 2 days after the 
burst ( |HaIpcrn et al. 2000| ). For this time scale to be consistent with an emission radius 
~ 10 16 cm, the relation t <; (1 + z)r /4cT 2 implies that the flow Lorentz factor must be <; 2. 
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However, given the comparatively high (~ 10 53 — 10 54 ergs) equivalent isotropic energy 
inferred for the emitting shock, it is unlikely that the Lorentz factor could become so low 
over such a relatively short distance (see bottom panels in Figs. f| and ^J). The problem is 
even more acute for GRB 970508, in which the X-ray emission feature detected 1 day 
after the burst again implies an emission radius ~ 10 16 cm (e.g., Lazzati et al. 1999), but 
where model fitting of the afterglow spectrum ~ 1 week after the GRB yields a radial 
scale ^ 3 x 10 17 cm along the line of sight to the center [see references in §[l|; this result 
is supported by radio scintillation measurements ([Frail ct al. 19971) 1. These values are 
mutually inconsistent, since the SNR shell could not have reached a distance of ^ 0.1 pc 
in one week even if it expanded at the speed of light. As we noted in §|3.1|, the afterglow 


emitting gas should decelerate rapidly after the forward shock encounters the SNR shell, 
and the shock transition into the shell would likely result in a discontinuous variation 
in the afterglow light curve. If the radius of the shell indeed coresponds to the value of 
R c j indicated by the X-ray emission-line data, then this is hard to reconcile with the fact 
that, in the case of GRB 970508, the light curve remained detectable and more or less 
smooth during a 450-day monitoring period, with the underlying flow evidently becoming 
nonrelativistic only after ~ 100 days ([Frail, Waxman, fc Kulkarni 2000| ; see also IGhcvalicil 

fc hi 200q) . 


The discrepancy between the inferred radius of the X-ray line-emitting shell and 
the deduced radial distance of the afterglow-emitting shock may be reconciled within 
the framework of the supranova/PWB scenario if the SNR shell and the PWB are not 


spherically symmetric. One possibility (suggested by Lazzati et al. 1999 and [Vietri et al. 


[2001| ) is that the supernova explosion does not eject matter along the SMNS rotation axis, 
where the GRB outflow is subsequently concentrated. An alternative possibility (which 
we discuss below) is that both the SNR and the PWB become elongated in the polar 
directions because of a preexisting density anisotropy in the GRB environment. Under these 
circumstances, a highly collimated GRB outflow (such as the one inferred in GRB 991216; 
Halpcrn ct al. 2000|) could reach a distance ^ 10 17 cm without encountering the SNR shell 
even as the lower-latitude regions of the shell (from which the X-ray line emission emanates) 
have radii <; 10 13 cm. In the case of GRB 991216, where the X-ray observations lasted 
between 37 and 40.4 hr after the burst (|Piro et ah 2000]) , the inferred effective spherical 
radius of the X-ray emitting shell (eq. |6l|) strictly corresponds only to angles 6 that lie 
in the narrow range ~ 27 — 28° (see eq. [58|]). If the jet opening half-angle is significantly 
smaller than these values and 74 is 10 16 cm at small values of 0, then the afterglow 
observations can in principle be consistent with the X-ray emission-line measurements.^ 
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The detection of an X-ray absorption feature would be compatible with this interpretation if it could be 
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The formation of a highly elongated PWB in the supranova scenario may be a natural 
outcome of the manner in which its environment was shaped by the progenitor star as 
well as of its intrinsic physical properties. The star that gave rise to an SMNS in a 
supranova event must have been massive, rapidly rotating, and magnetized, ft would have 
influenced the density distribution in its vicinity through episodes of strong mass loss, 
in particular during its red-supergiant and blue-supergiant evolutionary phases. There 
is strong observational evidence that the “slow” red-supergiant wind is often anisotropic 
(possibly as a result of fast-rotation and magnetic effects), transporting significantly more 
mass near the equatorial plane than in the polar regions. Subsequent stellar outflows that 
propagate into this mass distribution will assume an elongated morphology: this has been 
the basis of the “interacting stellar winds” class of models for the shapes of planetary 
nebulae (e.g., |Dwarkadas, Chevalier, fe Blondin 1996| ), in which the outflow represents the 
“fast” blue-supergiant wind, as well as of models of apparent SNR “protrusions,” in which 
the outflow corresponds to the supernova ejecta (e.g., [Blondin, hundqvist, fe Chevalier] 
1996). In these applications, the subsequent outflows have been taken to be effectively 


spherically symmetric. However, an even stronger collimation is achieved if the later 
outflow is itself anisotropic. In particular, if the fast wind is even weakly magnetized (with 
a dominant azimuthal field component), then, after passing through the wind shock where 
the field is amplified (an effect that will be especially strong if cooling is important behind 
the shock), the magnetic hoop stress will collimate the resulting interstellar bubble (e.g., 
Chevalier fc Luo 1994|) . In fact, as was argued by Gardiner & Frank (2001), the collimation 
may start even before the shock is encountered; this should be particularly pronounced 
in cases where magnetic stresses also play a dominant role in driving the fast wind (as in 
the Wolf-Rayet wind model of [Biermann fe Cassinelli 1993|) . The additional collimation 
provided by the magnetic field was suggested as the origin of strongly elongated planetary 
nebulae, which cannot be readily explained by purely hydrodynamic models. 

A pulsar wind expanding into the anisotropic density distribution created by the earlier 
(red-supergiant and blue-supergiant) stellar outflows will give rise to an elongated bubble 
(see, e.g., [Li fc Bcgclman 1992| for a discussion of PWB evolution in a stratified medium). 
Furthermore, since the pulsar wind is highly magnetized and cooling may be important 
in the supranova-induced PWB (see (| |2.3| ), the same magnetic collimation effects that are 
invoked in the modeling of planetary nebulae will act to increase the bubble elongation in 


demonstrated that the absorbing material was also located at a distance 10 16 cm from the irradiating- 
continuum source. In the only such case reported to date (GRB 990705), Lazzati et al. (2001) deduced a 
radius ~ 10 16 cm using a similar scheme to the one applied here to the interpretation of X-ray line emission. 
They have, however, also argued that the afterglow emission properties in this object may be consistent with 
a shock/SNR-shell encounter on this radial scale. 
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this case too. [In fact, the collimating effect of magnetic hoop stresses on plerionic supernova 
remnants was already discussed by Rees & Gunn (1974); it was subsequently modeled 
by Begelman & Li (1992).] Under these combined effects, it is quite plausible to expect 
that a bubble aspect ratio ^ 10 could be achieved, although this needs to be confirmed 
by an explicit calculation.^ Previous numerical simulations of outflows expanding into 
an anisotropic medium also make it likely that the column density of the swept-up SNR 
shell will be lower near the apex of the bubble than at larger values of 9, which should 
be relevant to the modeling of X-ray absorption and emission features as well as of the 
afterglow light curve. The expected departure of the PWB from sphericity might require 
a modification of the model presented in § |2.2| , which would probably be best done with 
guidance from numerical simulations. We nevertheless anticipate that the results obtained 
from the semianalytic model would remain at least qualitatively valid. Furthermore, if a 
strong elongation only occurs near the symmetry axis (which would be consistent with the 
data for GRB 991216 as well as with some of the existing numerical simulations), then 
even the quantitative predictions of the simple spherical model would still be approximately 
correct. 


4. Conclusions 

We propose to identify the environment into which afterglow-emitting shocks in at least 
some GRB sources propagate with pulsar-wind bubbles. Our results can be summarized as 
follows: 

• PWBs provide a natural resolution of the apparent difficulty of accounting for the 
high electron and magnetic energy fractions (e e and e^, respectively) inferred in a 
number of afterglow sources. This is because pulsar winds are expected to have a 
significant er component and to be highly magnetized. If high values of e e in fact 
prove to occur commonly in afterglow sources, then this would strengthen the case for 
a simple, “universal” explanation of this type. 

• An association of PWBs with GRBs is expected under several GRB formation 
scenarios, including the collapse of a massive star. In light of suggestive evidence that 
many of the afterglows observed to date may have a massive stellar progenitor, we 


33 In this case the width and centroid redshift of the observed X-ray emission lines may not be due entirely 
to the bulk motion of the SNR shell but may also have significant contributions from Compton broadening 


within the shell; see Vietri et al. 2001 
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have concentrated on this case. In particular, we considered the supranova scenario of 
VS, in which an intense pulsar-type wind from the GRB progenitor is a key ingredient 
of the hypothesized evolution. In this picture, the ejection of a highly energetic, 
ultrarelativistic pulsar wind is predicted to follow the supernova explosion and to last 
anywhere from several months to several years until the central object collapses to 
form a black hole, thereby triggering the burst. Recent detections of X-ray features 
in several GRB sources have been interpreted as providing strong support for this 
scenario. 

• To assess the implications of a PWB environment to afterglow sources in the context 
of the supranova scenario, we have constructed a simple, steady-state model of the 
early-time structure of a plerionic supernova remnant. We have been guided by 
Atoyan’s (1999) spectral modeling of the Crab, which yielded a lower initial wind 
Lorentz factor and a higher initial pulsar rotation rate than in previous estimates, 
and by other recent results on the Crab and Vela synchrotron nebulae, from which 
we inferred a plausible range of the wind magnetization parameter a w ( rs_/ 10“ 3 - 1). 
In contradistinction to previous models of the structure of plerionic SNRs, we 
have replaced the assumption that ideal MHD applies throughout the PWB with 
the postulate that the electromagnetic-to-thermal pressure ratio in the bubble 
remains constant after it increases to ~ 1. We have also explicitly incorporated 
synchrotron-radiation cooling. Although our solutions do not provide an exact 
representation of radiative (and thus intrinsically time-dependent) PWBs, we have 
verified that they generally do not depend on the detailed approximations that are 
adopted and are essentially characterized by a w and by a second parameter that 
measures the relative importance of radiative cooling within the bubble. It would be 
of interest to further develop this model and to investigate the possibility that it can 
be applied both to young radio pulsars and to GRB progenitors as members of the 
same general class of rapidly rotating and strongly magnetized neutron stars. 

• In view of the “hot” (relativistic) equation of state and high magnetization of the 
shocked wind, the effective hydrogen number density that determines the properties of 
a relativistic afterglow-emitting shock is given by nH, eq mv = [4 p+ (B 1 + £') 2 / An]/rn p c 2 , 
where B' and £' are, respectively, the comoving magnetic and electric fields and p is 
the particle pressure. For plausible values of the cooling parameter (and independent 
of the value of a w ), the derived values of rqpequiv span the density range inferred 
from spectral modeling of GRB afterglows. An interesting feature of the solutions is 
the predicted radial variation of nH,equiv within the bubble, which can mimic either 

a uniform-ISM or a stellar-wind environment, but which in general exhibits a more 
diverse behavior. Among other things, this model makes it possible to understand how 
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a GRB with a massive progenitor can produce an afterglow that shows no evidence of 
a stellar-wind or a high-density environment. 

• We have examined the dependence of the characteristic synchrotron spectral 
quantities in an afterglow-emitting shock that propagates inside a PWB on the 
bubble parameters and related them to the standard expressions derived under 
the assumption of a uniform-ISM environment. We found that, under typical 
circumstances, the standard expressions remain roughly applicable if one substitutes 
for e e , £b, and % their “equivalent” PWB expressions. We noted, however, that the 
parameter scaling laws would change in strongly radiative bubbles: these differences 
might be detectable in objects with high inferred ambient densities. 

• Finally, we considered the possible observational manifestations of the dense supranova 
shell that surrounds the PWB in this picture. In particular, we discussed how the 
X-ray emission features detected in objects like GRB 991216 may be interpreted 

in the context of a supranova-generated PWB. We concluded that both the X-ray 
features and the afterglow emission could be explained by this model if the PWB were 
elongated, and we argued that such a shape might be brought about by anisotropic 
mass outflows from the GRB progenitor star. 
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Fig. 1.— The wind-shock radius R s (top panels) and the equipartition radius R cq (bottom 
panels), normalized by the bubble radius Rb, as functions of the parameter cq (defined in eq. 
HH) for 4 values of the pulsar-wind magnetization parameter a w . The left panels present 
solutions obtained under the particle-conservation constraint (eq. ||4|]), whereas the right 
panels show solutions derived by fixing the gas speed at Rb (eq. [j45fl ). 
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Fig. 2.— The panels show the radial distributions of /3 (the gas speed in units of the speed of 
light), n (the particle number density, for which the dimensional scaling corresponds to the 
fiducial values of the model parameters), y e (the random electron/positron Lorentz factor, 
normalized by the pulsar-wind Lorentz factor j w ), p (the particle pressure, normalized by 
Pi ; see eq. |37j), 0 = (1 + S' / B') 2 B' 2 /Sirp, and S' / B' (the ratio of the comoving electric 
and magnetic fields) in the model PWBs for 4 values of the pulsar-wind magnetization 
parameter a w (listed at the top of the respective columns of panels). The solution curves in 
each panel correspond to 5 values of the parameter ap. 0.45, 5.04, 45.3, 504, and 4539 (for 
endpoints progressing respectively from left to right); they were obtained by imposing the 
particle-conservation constraint (eq. 
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Fig. 3.— Same as Fig. |2]. except that the solutions were obtained by using the final-speed 
constraint (eq. [[RJ) instead of the particle-conservation constraint. 
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Fig. 4.— The panels show the radial distributions of n H ,equiv (the effective hydrogen number 
density, eq. [|TJ), 2 — k (with k = — dlog nH,equiv/dlogr being the effective power-law index of 
the equivalent hydrogen density distribution), e eieq uiv (the equivalent electron energy fraction, 
eq. [^2|), e_B,equiv (the equivalent magnetic energy fraction, eq. |53[), incorrect = ^-/'^H.cquiv (the 
flux correction factor, eq. |54|]), and the Lorentz factor of the shocked bubble material (for a 
spherical shock driven into the PWB by an outflowing mass of energy E and initial Lorentz 
factor r 0 ) for the PWB solutions depicted in Fig. ||. The displayed results correspond to 
the fiducial values of the model parameters. The solid and dash-dotted curves in the bottom 
panels correspond to E = 10 52 and 10 53 ergs, respectively. 
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Fig. 5.— Same as Fig. [|, except that the results correspond to the PWB solutions depicted 
in Fig. [|. 































































